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Abstract 



We present a search for neutrino-induced events containing a single, exclusive pho- 
ton using data from the NOMAD experiment at the CERN SPS where the average 
energy of the neutrino flux is ~ 25 GeV. The search is motivated by an excess of 
electron-like events in the 200-475 MeV energy region as reported by the Mini- 
BOONE experiment. In NOMAD, photons are identified via their conversion to 
e+e" in an active target embedded in a magnetic field. The background to the sin- 
gle photon signal is dominated by the asymmetric decay of neutral pions produced 
either in a coherent neutrino-nucleus interaction, or in a neutrino-nucleon neutral 
current deep inelastic scattering, or in an interaction occurring outside the fiducial 
volume. All three backgrounds are determined in situ using control data samples 
prior to opening the 'signal-box'. In the signal region, we observe 155 events with 
a predicted background of 129.2 ± 8.5 ± 3.3. We interpret this as null evidence 
for excess of single photon events, and set a limit. Assuming that the hypothetical 
single photon has a momentum distribution similar to that of a photon from the 
coherent vr'' decay, the measurement yields an upper limit on single photon events, 
< 4.0 X lO"'^ per i/^ charged current event. Narrowing the search to events where 
the photon is approximately collinear with the incident neutrino, we observe 78 
events with a predicted background of 76.6 it 4.9 it 1.9 yielding a more stringent 
upper limit, < 1.6 X 10~^ per charged current event. 

Key words: single-photon neutrino neutral current coherent pion 
PACS: 13.15.+g, 13.85.Lg, 14.60.Lm 



1 The Question 



The MiniBooNE experiment, in the neutrino-mode, has reported an excess of 
129±43 (Stat ©Syst) electron-like events in the 200-475 MeV energy range [1], 
a range largely unaffected by the t/^ — >■ Ue oscillations at Am^ > 1 eV^, which 
is expected to dominate the 475-1250 MeV range [2]. Assuming that the ef- 
ficiency of these events is 0.25, similar to that of the i/g-induced electrons, 
the rate of the excess with respect to i^^-CC is ~ 3 x 10"^. In the antineu- 
trino mode, initially the MiniBOONE data were consistent with background in 
the 200-475 MeV range [3] . Recently, with additional antineutrino data, they 
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report a 2a excess in this region [4]. MiniBOONE lacks the resohition to deter- 
mine if the excess is due to e^, e+, or a photon. In neutrino interactions single 
photons do not occur: they occur in pairs from the 7r° 77 decay. Evidence 
of events with a single photon and nothing else in neutrino induced neutral 
current (NC) interaction will signal either new or unconventional physics. This 
anomaly needs to be checked experimentally. 

The low-energy MiniBOONE excess has motivated novel hypotheses invoking 
conventional, to unconventional, to entirely new physics processes. An expla- 
nation [5] of the excess attributed to the internal bremsstrahlung associated 
with muons in i>'^-induced charged current interaction (CC) has been refuted 
by MiniBooNE [6]. Harvey, Hill, and Hill [7,8] postulate a new anomaly- 
mediated interaction between the //-induced Z°-boson, the photon, and a 
vector meson, such as a;, coupling strongly to the target nucleus. A similar, 
single photon process has been envisioned by Jenkins & Goldman [9]. In the 
anomaly-mediated neutrino-photon (ANP) interaction, vM — >■ vM^, where M 
is the target nucleus or nucleon, the observable is a single 7 , approximately 
collinear with the incident neutrino direction in the lab frame with a recoil- 
ing A/" which remains intact and largely undetected; the form factor of J\f is 
expected to induce the single 7 to be pulled forward. The ANP cross-section 
is expected to increase with neutrino energy {E^,). The ANP interaction is 
akin to the exclusive tt^ production when the neutrino coherently scatters off 
a target nucleus (Coh7r° ), -|- M v -\- Af -\- 7r°, where the only observ- 
able is 77° — >■ 77. The phenomenology of ANP developed by Hill in [10], and 
largely focussed on the low energy neutrino interaction {E,y ~ GeV), 
predicts a cross-section sufficiently large to explain the MiniBOONE anomaly. 
If true, the ANP interaction will also impact astrophysical processes such as 
neutron star cooling. Hypotheses involving new physics include: decay of a 
heavy neutrino into a light neutrino and a photon, i^/j — >• 1/ -|- 7, by Gni- 
nenko [11,12]; CP- violation in 30 2 model by Maltoni & Schwetz [13], and by 
Goldman et al. [14]; extra-dimensions in 301 model by Pas et al. [15]; and 
models by Giunti and Laveder [16]; Lorentz- violation by Katori et al. [17]; 
CPT-violation in 301 model by Barger et al. [18]; new gauge bosons with 
sterile neutrinos by Nelson & Walsh [19]; and soft-decoherence by Farzan et 
al. [20]; etc. The NOMAD collaboration has reported a search of heavy-iv 
which mixes with i/^-'s produced in the SPS proton-target, and then decays 
into e~e^ pair in the detector [21]. The search focussed on highly collinear 
(C = [1 - cosGee] < 2 X 10"^) and high energy {E^^ > 4 GeV) e-e+ 
events, where Q^e and E^e are the polar angle and the energy of the con- 
verted photon, and resulted in a single event consistent with the estimated 
background. The present work extends this search with no cut on C. 

The high resolution NOMAD data allow a sensitive search for neutrino in- 
duced single photon events (I7 ). The detector measures the energy and emis- 
sion angle of the photon. Additionally, the detector affords the redundancy to 
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measure in situ all relevant backgrounds relieving the reliance on Monte Carlo 
simulation (MC) of the conventional processes. As regards the MiniBooNE 
low-energy anomaly, we note that the average energy of the SPS neutrino 
flux at NOMAD, ~ 25 GeV, is much higher than that of MiniBOONE, 
E,y ^ 1 GeV. However, the NOMAD data provide a stringent check in a dif- 
ferent energy range. If a mechanism explaining the anomaly were applicable 
to higher energies, such as ANP or ly-decay hypotheses, then this search, even 
if negative, will furnish meaningful limits. 



2 Beam and Detector 

The Neutrino Oscillation MAgnetic Detector (NOMAD) experiment at CERN 
used a neutrino beam produced by the 450 GeV protons from the Super Proton 
Synchrotron (SPS) incident on a beryllium target and producing secondary 
TT^, K^, and mesons. The positively charged mesons were focussed by two 
magnetic horns into a 290 m long evacuated decay pipe. Decays of tt^, K^, 
and produced the SPS neutrino beam. The average neutrino flight path to 
NOMAD was 628 m, the detector being 836 m downstream of the Be-target. 
The SPS beamhne and the neutrino flux incident at NOMAD are described 
in [22]. The i/-flux in NOMAD is constrained by the tt^ and produc- 
tion measurements in proton-Be collision by the SPY experiment [23,24,25] 
and by an earlier measurement conducted by Atherton et al. [26]. The Ei,- 
integrated relative composition oi Vf^.V f^:Ue.Ve CC events, constrained in situ 
by the measurement of CC-interactions of each of the neutrino species, is 
1.00 : 0.025 : 0.015 : 0.0015. Thus, 97.5% of the events are induced by 
neutrinos with a small anti-neutrino contamination, similar to that reported 
by MiniBOONE [2]. 

The NOMAD apparatus, described in [27], was composed of several sub- 
detectors. The active target comprised 132 planes of 3 X 3 m^ drift chambers 
(DC) with an average density similar to that of liquid hydrogen (0.1 gm/cm^). 
On average, the equivalent material in the DC encountered by particles pro- 
duced in a iv*- interaction was about half of a radiation length and a quarter of a 
hadronic interaction length (A). The fiducial mass of the NOMAD DC-target, 
2.7 tons, was composed primarily of carbon (64%), oxygen (22%), nitrogen 
(6%), and hydrogen (5%) yielding an effective atomic number, A. =12.8, simi- 
lar to carbon. Downstream of the DC, there were nine modules of transition ra- 
diation detectors (TRD), followed by a preshower (PRS) and a lead-glass elec- 
tromagnetic calorimeter (ECAL). The ensemble of DC, TRD, and PRS/ECAL 
was placed within a dipole magnet providing a 0.4 T magnetic field orthog- 
onal to the neutrino beam line. Two planes of scintillation counters, Ti and 
T2, positioned upstream and downstream of the TRD, provided the trigger 
in combination with an anti-coincidence signal, V, from the veto counter up- 
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stream and outside the magnet. Downstream of the magnet was a hadron 
calorimeter, followed by two muon-stations each comprising large area drift 
chambers and separated by an iron filter. The two stations, placed at 8- and 
13-A's downstream of the ECAL, provided a clean identification of the muons. 
The schematic of the DC-tracker with a 7-conversion in the Y-Z view is shown 
in Figure 1. The charged tracks in the DC were measured with an approximate 
momentum (p) resolution of CTp/p = 0.05/\/L 0.008p/\/L^ {p in GeV/c 
and L in meters) with unambiguous charge separation in the energy range of 
interest. The experiment recorded over 1.7 million neutrino interactions in the 
active drift-chamber (DC) target in the range 0{1) < < 300 GeV. 
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Fig. 1. Schematic of the DC tracker with a single 7-conversion event in NOMAD 
in Y (vertical) vs Z (horizontal) view. The energy and angle with respect to the 
neutrino direction of the 7 are 3.1 GeV and 0.095 Radians. 

3 The Exclusive Photon Signature 

In NOMAD the cleanest signature of a 7 is a V° arising from its conversion 
in the DC into e~ and 6+ tracks as shown in Figure 1. The 7 -momentum is 
reconstructed by measuring the 3-momenta of the associated e~ /e'^ tracks in 
the DC. The ECAL measures energy associated with the eT je^ tracks plus 
any unconverted photons. Because the trigger counters are situated at the 
downstream end of the detector, only a fraction of events containing single 
photon will trigger the apparatus, namely events where the photon undergoes 
conversion in DC and the resulting e+ or eT traverses the trigger counters. 
Events with unconverted photons reaching the ECAL will not trigger the 
apparatus. For example, about 29% of the CohTr'' events, containing a tt" 
and 'nothing else', trigger the apparatus; the loss arises from the unconverted 
photons and, among the converted photons, from the e~/e"'" tracks that do 
not reach the downstream trigger counters. 

Single photon events will manifest as a V° , composed of e~ and 6+, with 
no additional energy in the ECAL. A converted photon is defined as a V° 
such that the invariant mass of the e~ and 6+ [M^e) is less than 100 MeV 
which selects the converted photons with 95% purity and 97% efficiency. Fur- 
thermore, if the single photon is approximately coUinear with the incident 
neutrino then the coUinearity of the photon-track offers an additional kine- 
matic handle. Because we do not have a definitive simulation for the single 
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7 process, we use one of the photons from the Cohn^ , while ignoring the 
other photon, as a guide for the signaL We define two signal boxes. Boxl 
comprises single •y events with negligible energy in the ECAL associated 
with neutral particles ( 7 , neutron, etc.) quantified by the quantity, PAN = 
[E^ — ENeut] I [-E'T + -E'iVeut], > 0.9, where E^ is the energy of the photon 
derived from the e~ /e'^ tracks and En^ui is the ECAL energy associated 
with other neutral particles, i.e. not associated with / e~ tracks. (For cal- 
culation of Ejsfeut see, for example, [28].) The PAN > 0.9 cut is model 
independent. Box2, a subset of Boxl, comprises single 7 events provided 
that they are approximately coUinear with the incident neutrino, character- 
ized by the variable, C — -^7 [1 ~ cos{0^)] < 0.05. where E^ and 6^ are 
the photon's energy and polar angle with respect to the neutrino beam. (The 
variable C has the property that its distribution depends weakly on the inci- 
dent neutrino energy.) Approximately 90% of the photons from CohTT*' have 
C < 0.05. Although the C ^ 0.05 cut will be somewhat model-dependent, at 
the SPS energies {E^, ~ 25 GeV) most exotic models suggest that the single 
photon will be approximately coUinear with the incident neutrino. 

The backgrounds to the single 7 are all interactions where all daughter 
particles, except one photon, evade detection. The background includes the 
CohTT^ interaction, the deep inelastic tz-NC scattering (NC-DIS), and u- 
interactions occurring outside the fiducial volume (OBG). Note, that with a 
single converted 7 in DC the coordinates of the i/-interaction cannot be 
ascertained. The sample with PAN < 0.9 is composed of events with a V° 
and neutral energy in the ECAL as expected from a 7r° decay and constitutes 
an important control sample. All three backgrounds are determined in situ as 
described below. 

The analysis proceeds as follows. First, the sample with a single photon- 
conversion in DC is selected. Second, the three backgrounds to the single 
photon sample are calibrated using data outside Boxl. The analysis leads 
to an experimentally constrained prediction of the background events in the 
signal region, Boxl and Box2. Finally, boxes are opened and the predictions 
are compared with the observed data. 



4 Selection of Single Photon Events 

We select t'-induced events where a single photon converts (V° ) within the 
fiducial volume of the DC target. The analysis uses the entire NOMAD data 
and the associated Monte Carlo (MC) samples described in [29]. The NC- 
DIS sample, defined by requiring that the generated invariant hadronic mass 
squared (W^) be > 1.96 GeV^, is normalized to 0.53 X 10^ events which 
corresponds to 37% of the i^j^-CC. The normalization of the NC-Resonance 
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3,076 


10,547 


Eee > 1.5 GeV 
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69 


1,846 


4,543 



Table 1 



Preselection of events with a single, converted 7 in Colivr , NC-Resonance, and 
NC-DIS MC samples, and data. The MC samples, having a much larger statis- 
tics than data, are normalized to the expectation as described in the text. The 
|(X, Y — 5)|pROJ refer to X- and Y-position of the 7 -vector when projected to 
the most upstream DC (Zmiii)- 

(W^ < 1.96) sample is set at 3.5% of the NC-DIS. The CohTr" interaction is 
simulated using the Rein-Sehgal (RS) model [30] which agrees well with our 
measurement [31]. 

The MC and data samples are subjected to a preselection requiring: (a) the 
presence of a converted photon whose reconstructed conversion point (X, Y, 
Z) be within the fiducial volume, |X, (Y — 5)| < 130 cm and 
ZMin < Z < 405 cm where Zmiii is 5 and 35 cm for the two configurations 
of the detector composing more than 95% of the NOMAD data; (b) no recon- 
structed muon (/ilD-veto) and a single photon-induced V" ; (c) the X and Y 
coordinate of the 7 -vector when extrapolated back to the upstream most DC 
at ZMin be within the transverse fiducial volume, |X, Y — 5|proj < 130cm, 
largely eliminating 7's that enter from the sides; and (d) the energy of the 
e~e+ pair be > 1.5 GeV which reduces the NC-DIS and OBG background 
while having a small effect on the CohTr" sample which also serves as a guide 
for the single 7 signal. The preselection reduces the NC-DIS and data sam- 
ples by more than a factor of one hundred. It is noteworthy that only a small 
fraction of the NC-Res pass the selection since most of the photons, coming 
from 7r° emitted at large angles, either fail to trigger the apparatus or have 
energies below 1.5 GeV. Given the paucity of the NC-Resonance events, it 
has been added to the NC-DIS sample. Finally a negligible fraction of Uf^-CC 
(< 10~^) pass the selection; consequently the CC sample has been ignored in 
this analysis. The preselection is presented in Table 1. 

The final event selection follows the preselection with more stringent require- 
ments. The vertex coordinates of V° are required to be within 
|X, (Y — 5) I < 120 cm and ZMin ^ Z < 405 cm. Two additional cuts (Clean- 
V° ) are imposed to reduce outside background by requiring that there be no 
tracks upstream of the photon conversion and that there be no hits associ- 
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Cut 


CohTrO -RS 


NC-DIS © Resonance 


Start 


516 


1915 


Tighter Fid-Cuts 


483 


1775 


Mee < 50 MeV & Clean-V° 


386 


400 



Table 2 

Final selection of single 7 Events in the MC Samples: The MC samples have been 
normalized as presented in Section 4. 

ated with the tracks composing the converted- 7 in the most upstream DC. 
Finally, the Mee cut is tightened to < 50 MeV which increases the photon 
conversion purity to > 98% while reducing the efficiency to 93%. Table 2 
summarizes the selection of events in the MC samples. The preselected data 
are subjected to identical cuts. Distributions of the X, Y, and Z coordinates 
of the I7 vertex in Data and the corresponding MC prediction, composed of 
Coh7r° OBG, NC-DIS, agree. Only the control sample with PAN < 0.9 is 
examined to check and constrain the background prediction. The calibration 
of the three backgrounds is presented next. 



5 Constraining the Backgrounds and the PAN < 0.9 Sample 

The prediction of the backgrounds to the single 7 sample is data driven. The 
estimation of the CohTr'^ induced I7 needs to be based on the observed 2V° 
sample where both photons convert in the DC. Furthermore, Monte Carlo 
simulations can neither reliably provide the OBG nor the NC-DIS contribu- 
tion to the I7 samle. These backgrounds need to be determined using data 
themselves. 

First, we present the calibration of the Coh7r° contribution to the I7 sample 
because it is the simplest. The analysis of the 2V° sample, where two photons 
and nothing else is observed, provides a measurement of the Coh7r° process. It 
yields a rate with respect to the inclusive i>'^-CC, [3.21 ± 0.46{stat syst)] X 
10~^ [31]. The measured kinematic distributions of the 27 data are consistent 
with the RS-Coh7r" model. Consequently we use the simulated RS-Coh7r° 
events, normalize it to the measured cross-section, subject these to the I7 
analysis, and obtain the Coh7r° contribution to the I7 sample. The analysis 
yields a normalized prediction of Cohiv^ -induced I7 with a 14.5% precision, 
shown in Table 4. 

Next, we present the calibration of the OBG contribution to the I7 sample. 
This, too, is accomplished using the 2V° sample as in [31]. The two recon- 
structed photon momentum vectors enable one to determine the tz-interaction 
vertex by extrapolating the vectors upstream and finding the coordinates of 
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DCA-Z>Zmin 


DCA-Z<Zmm 


27 Data 
27 OBG-Data 


381 
451 


169 
927 



Table 3 



Events passing and failing the Z-cut in 27 samples in Data and OBG-Data 

their distance of closest approach (DCA). The procedure defines the DCA- 
vertex with coordinates denoted as DCA-X, DCA-Y, and DCA-Z. The DCA- 
vertex resolution is well understood. Among the 27 sample 169 events, out of 
550 events (see Table 3), have a DCA vertex of origin upstream of the ZMin 
providing a calibration for the OBG-I7 as follows. A different data sam- 
ple, OBG-Data, which includes charged tracks, is selected in which the ver- 
tex reconstructed using these charged tracks is upstream of the fiducial limit 
(Z < ZMin)- In this sample, the charged tracks are ignored and the events 
subjected to the 27 analysis (see [31]). This OBG-Data sample with 27 
yields 927 events with DCA-Z< Zmiu and 451 events with DCA-Z> Zmiii 
as shown in Table 3. The number of 27 events with DCA-Z> ZMinj orig- 
inating from events occurring at Z < ZMin but producing no visible tracks, 
is then calculated as (451/927) X 169 = 82.2 ± 6.1. This constitutes the 
OBG background to the 27 sample. The so normalized OBG sample is then 
subjected to the I7 analysis yielding a calibrated OBG prediction, with a 
7.7% error (driven by the statistics of 169 events), shown in Table 4. 

Finally, we present the cahbration of the NC-DIS contribution to the I7 sam- 
ple using the control region PAN < 0.9 dominated by events containing a 
photon-conversion accompanied by a neutral energy cluster in the ECAL. Fig- 
ure 2 shows the observed PAN distribution for the full I7 sample although 
the signal region PAN > 0.9 is not looked at till all backgrounds are final- 
ized. The figure evinces agreement between data and the prediction for PAN 
< 0.9. Furthermore, the observed kinematic variables associated with the 
neutral cluster in ECAL are found to be consistent with the prediction. The 
observed momentum distribution of the photon is consistent with the predic- 
tion as evidenced by Figure 3. Lastly, Figure 4 shows that the measured C-y 
distribution agrees with the MC prediction. Because the tt" -induced photons 
in NC-DIS will have a broader ^-y than those in Coh7r° , we use Ct> 0.05 
to normalize the NC-DIS. The NC-DIS normalization factor is 1.13 ih 0.14. 
(This normalization factor has been applied to the NC-DIS in the figures.) 
The logic behind using this control sample to constrain the background in 
the signal region is that I7 events come from 7r"-dominated neutral current 
events where one of the photons evades detection. One concern, however, is 
that of a possible systematic bias when using the PAN < 0.9, a region where 
both photons are reconstructed in DC and ECAL, i.e. the two photons are 
emitted in the forward direction, to predict the NC-DIS in the PAN > 0.9 
where the second photon misses the ECAL. To check this concern, ordinary 
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NC events, where charged tracks define the event, and, thereby, reheving the 
forward bias in the I7 reconstruction due to the trigger requirement, are 
selected in data and MC. A single I7 is selected downstream of the event. 
The charged tracks are ignored, the I7 is subjected to the current analysis. 
The Data/MC ratio is examined in and C-y variables. The check reveals 
that the shapes of these variables in MC are consistent with those of data: the 
ratio is unity within ±5% in the kinematic range, well covered by the 12% 
error ascribed to the NC-DIS. Table 4 presents the NC-DIS contribution to 
the I7 sample. 

To sum up, we have a calibrated prediction of Coh7r° , OBG, and NC-DIS 
contributions to the I7 sample. Table 4 lists the corresponding errors, all 
statistical in nature, as determined by the respective data control samples. 
Additionally, we estimate the error in the I7 reconstruction to be < 2.5%, 
which has a negligible effect on this analysis. As a final check, we reproduced 
the results presented in the [21]. This exercise illustrates that the current 
analysis is two orders of magnitude less stringent than that in [21] which was 
focussed on the search for a heavy neutrino coupling to u^. 
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Fig. 2. Comparison of PAN , the energy asymmetry between the I7 -momentum 
and the neutral energy in the ECAL, between Data (Symbol) and prediction: CohTr'' 
(Blue-hatched), OBG (Green-dash-dot), NC-DIS (Red-dash), and total (Black-his- 
togram). The signal region, PAN > 0.9, is not looked at till the analysis is com- 
plete. 
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Fig. 3. Comparison of -Ry between data and MC in PAN < 0.9 region. 
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Fig. 4. Comparison of C-y distribution of the I7 sample with PAN < 0.9 in data 
and MC. The C-y > 0.05 region is used to normahze the NC-DIS. 

6 Results 



Figure 1 shows a representative event in the signal region, PAN > 0.9. 
The first signal-box, Boxl, with PAN > 0.9, yields 155 observed events 
with a predicted background 129.2 ih 8.5 ih 3.3 events yielding an excess of 
25.8 ih 15.5 events. Figure 5 and Figure 6 present the Pj and (^j comparison 
between data and MC. Absence of significant excess leads us to interpret this 
as null evidence for excess of single photon events. Assuming the error on 
the background-prediction to be Gaussian, we derive an upper limit of < 51 
events at 90% CL in Boxl. 

The second signal-box, Box2, with PAN > 0.9 and 0.05, yields 78 ob- 
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Cut 


Coh7r° -RS 


NC-DISffiRes 


OBG 


Total 


Data 


V odlllUiC 






341 3 


1 1 97 3 


1 1 zLQ 


TV l\ Tj'-MMrtW 

iviC/ H/rror 






7 7% 
< . 1 /o 






Background 












PAN ^ n q 


353.9 


317.7 


296.5 


998 1 +69 9+'^5 f) 


994 


Signal 












PAN > 0.9 


32.0 


52.4 


44.8 


129.2±8.5±3.3 


155 


PAN > 0.9 e C7< 0.05 


22.8 


22.6 


31.2 


76.6±4.9±1.9 


78 



Table 4 



The I7 Sample: Presented are the normalized Coh7r° , NC-DIS, and OBG predic- 
tions for the I7 sample. The MC errors of the three components are all statistical 
in nature, as determined by the respective control samples. The systematic error 
due the V° -reconstruction is shown under the 'Total' column. Data are shown in 
the last column. 

served events with a predicted background 76.6 ih 4.9 ±1.9 events yielding 
an excess of 1.4 it 10.3 events. Events in the Box2 exhibit kinematic charac- 
teristics consistent with the background prediction. The vertex-distributions 
of the I7 events in the signal-region are in agreement with the MC predic- 
tion; within errors, the and C,-^ distributions, shown in Figure 7 
and Figure 6. In addition, the observed coUinearity (C) of the photon matches 
that of the prediction, see Figure 8. Assuming the error to be Gaussian, we 
derive an upper hmit of < 18 events at 90% CL. Table 4 presents the final 
enumeration. 



The only remaining task is to set an upper limit on the rate of single 7 events. 
To determine the signal efficiency, we assume that the 'signal' photon has 
kinematics similar to that of one of the photons from the CohTr^ interaction 
(the other photon is removed from the simulation). Using the RS Coh7r° - 
model, the signal efficiency for Boxl, PAN > 0.9, is 8.8%. The efficiency for 
Box2, PAN > 0.9 and C < 0.05, is 8.0%. The number of fully corrected 
v^-QQ in the same fiducial volume is measured to be 1.44 X 10^. We obtain 
the following upper limits on the rate of single photon events in ^'-interactions: 



o-(Single - 7) ^ , ^ _ .,„_4 



< 4.0 X 10-^(90%CL) (1) 



< 1.6 X 10-^(90%CL) (2) 



cr (Single, Forward — 7) ^ ^ ^ -4 
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In summary, we have presented a search for single photon events in interactions 
of neutrinos with average energy ~ 25 GeV. All relevant backgrounds are 
constrained using data control samples. No significant excess is seen. Assuming 
that the hypothetical signal has kinematics similar to those of a photon from 
the CohTT^ interaction, the analysis imposes an upper limit on the rate of 
excess of single photon events of < 4.0 x 10~^ per u^-CC at 90% CL; with 
an additional soft collinearity cut (permitting 90% of 7 from Coh7r° ) the limit 
is < 1.6 X 10~^ per u^-CC at 90% CL. Following the report on superluminal 
neutrinos by the OPERA collaboration [32] , we are conducting a speciahzed 
search for very forward e~e'^ pairs. 
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Fig. 5. Comparison of P-y between data and MC in Boxl, PAN > 0.9 region. 
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Fig. 6. Comparison of C-y distribution between data and MC in Boxl, PAN > 0.9 
region. 
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Fig. 7. The Momentum of I7 in Box2: The observed distribution is consistent with 
the MC-prediction. 
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Fig. 8. Consistent collinearity of ly in data and MC in the Box2. 
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